Superstatistical distributions from a maximum entropy principle 
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We deal with a generalized statistical description of nonequilibrium complex systems based on 
least biased distributions given some prior information. A maximum entropy principle is introduced 
that allows for the determination of the distribution of the fluctuating intensive parameter j3 of a 
superstatistical system, given certain constraints on the complex system under consideration. We 
apply the theory to three examples: The superstatistical quantum mechanical harmonic oscillator, 
the superstatistical classical ideal gas, and velocity time series as measured in a turbulent Taylor- 
Couette flow. 
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I. INTRODUCTION 

Many complex systems in physics, biology, medicine, 
and economics exhibit a spatio-temporally inhomoge- 
neous dynamics that can be effectively described by a su- 
perposition of several statistics on different time scales, in 
short a 'superstatistics' % S S i, 1J, 0, & 1, M M 03 • 
The concept of such a superposition of statistics was first 
systematically discussed in [![, in the mean time many 
applications for a variety of complex systems have been 
pointed out [H Q 03 M, 03, OS 01 [H . Essential for 
this approach is the existence of an intensive parameter (3 
that fluctuates on a much larger time scale than the typ- 
ical relaxation time of the local dynamics. In a thermo- 
dynamic setting, f3 can be interpreted as a local inverse 
temperature of the system, but much broader interpre- 
tations are possible. Locally, the system is described by 
equilibrium statistical mechanics with inverse tempera- 
ture f3, whereas globally there is yet another statistics 
of the inverse temperature (3. The two effects produce 
a superposition of two statistics, or in a short, a 'su- 
perstatistics'. Related statistical tools play of course an 
important role in the theory of stochastic processes, see 
e.g. 011, 11, [11. 

The stationary distributions of superstatistical sys- 
tems, obtained by averaging over all /3, typically exhibit 
non-Gaussian behavior with fat tails, which can decay 
e.g. with a power law, or as a stretched exponential, or in 
a more complicated way In general, the superstatis- 
tical parameter (3 need not to be an inverse temperature 
but can also be interpreted as an effective friction con- 
stant, a changing mass parameter, a changing amplitude 
of Gaussian white noise, the fluctuating energy dissipa- 
tion in turbulent flows, a fluctuating volatility in finance, 
an environmental parameter for biological systems, or 
simply a local variance parameter extracted from a given 
experimental time series. Recent applications of the con- 
cept include hydrodynamic turbulence I2I 12f3. 125ll2t| . pat- 
tern forming systems [l3[, cosmic rays solar flares 
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(TH , share price fluctuations random ma- 

trix theory [l?], HI , random networks |3l| , multiplicative- 
noise stochastic processes |32| . quantum systems at low 
temperatures @, wind velocity fluctuations [HI], hydro- 
climatic fluctuations [III] , the statistics of train departure 
delays [33| and models of the metastatic cascade in can- 
cerous systems [34] . 

In equilibrium statistical mechanics it is clear how to 
obtain the relevant probability distributions describing 
the long-term behavior of the system under considera- 
tion. These are the canonical distributions and they fol- 
low from a maximum entropy principle. However, su- 
perstatistical systems are nonequilibrium systems with a 
stationary state which is a mixture of canonical distribu- 
tions. It is a priori not clear how to obtain the mixing 
distribution of the fluctuating parameter from first prin- 
ciples. A promising idea to tackle this problem is to de- 
velop a more general type of thermodynamics for super- 
statistical systems which leads to a generalized maximum 
entropy principle that fixes these distributions. Early at- 
tempts in this direction were made by Tsallis and Souza 
^1 and later by Abe et al [Hj], Crooks [3(| and Naudts 
37]. Inspired by these early considerations, in this paper 
we develop a generalized formalism that is a) conceptu- 
ally simple b) applicable to both, classical and quantum 
systems c) consistent with experimental observations. As 
a result, we obtain a statistical theory that can be ap- 
plied to a large variety of complex systems and which 
further develops the earlier ideas of Abe, Beck, Cohen, 
Crooks and Naudts. 

This paper is organized as follows. In section 2 we 
clarify our notation and recall the basic concept of time 
scale separation that lies at the heart of any superstatisti- 
cal description. In section 3 we introduce our generalized 
maximum entropy principle and discuss the relation be- 
tween our formalism and the previous approaches of Abe, 
Beck, Cohen, Crooks and Naudts. In section 4 we dis- 
cuss some physically relevant conditions on the relevant 
class of probability densities. In the following sections 
we apply our theory to three important examples: The 
superstatistical quantum mechanical harmonic oscillator 
(section 5), the superstatistical ideal gas (section 6) and 
velocity fluctuations as observed in a turbulent time se- 
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ries (section 7). Our concluding remarks are given in 
section 8. 



II. BASIC CONCEPTS 

The crucial assumption of superstatistics is that 
the statistical description of certain classes of complex 
nonequilibrium systems can be split into two levels that 
have a large time scale separation. The total system is 
divided into spatial cells, each in local equilibrium, but 
the temperatures of the different cells don't have to be 
equal. As a consequence, in very good approximation the 
local properties of the different cells can be described us- 
ing the standard Boltzmann-Gibbs formalism. The main 
problem is then the determination of the distribution of 
the temperature at the higher level of the total nonequi- 
librium system. Clearly, the Boltzmann-Gibbs formalism 
is not applicable at this level. 

Locally, in each cell the average (A)h of an observ- 
able A is calculated with respect to the Boltzmann-Gibbs 
probability measure 

where (3 is the inverse temperature, H is the Hamiltonian 
that describes the properties of each spatial cell of the 
system, and Z((3) is the partition function. In classical 
statistical mechanics, p(H ; (3) is a probability distribu- 
tion and the local average (A)h is defined by 

(A) H = J dTp(H;(3)A, (2) 

with r being the phase space. In quantum statistical 
mechanics, p(H ; (3) is a density operator and the local 
average {A)h is defined by 

(A) H = Trp(H;/3)A, (3) 

with H and A being operators acting on the correspond- 
ing Hilbert space. We introduce the following shorthand 
notation for the local energy E((3) and local entropy S((3) 

E(f3) = (H) H and S(p) = - (In p(H ; f3))„. (4) 

From a thermodynamic point of view, the Hamiltonian is 
an observable and the temperature is the corresponding 
control parameter (intensive variable). By measurement 
of the average value of the observable one can estimate 
the value of the corresponding control parameter. We are 
interested in the statistical average of an observable A 
of the total nonequilibrium system which has a different 
inverse temperature in each cell. For this global average 
we will use following notation 

/>OC 

((A) H ) fi = / d(3f{(3;\i)(A) H . (5) 
Jo 



Here /(/?; Aj) is the probability density of (3 in the var- 
ious spatial cells, which depends on a set of parameters 
{Xi} (in our notation we suppress the brackets {}). The 
parameters A^ can be interpreted as the control parame- 
ters corresponding with some measurable nonequilibrium 
observables. Our goal in the following is to find a general 
principle for the determination of f((3; Xi), given certain 
information that we have on the complex system. 

III. MAXIMUM ENTROPY 

Let us first recall the maximum entropy principle for 
equilibrium statistical mechanics, after that we will pro- 
ceed to the superstatistical extension. An impressive 
amount of experimental results shows that assuming the 
Boltzmann-Gibbs distribution as the equilibrium distri- 
bution of a system is a very good approximation. In- 
formation theory gives a deeper understanding to this 
success [38] . Usually, the only experimental information 
that is available about a system is the average value of 
some observables. Therefore, it is natural to use the least 
biased distribution, given this prior information, as the 
equilibrium distribution of the system. The practical tool 
to obtain this least biased distribution is the maximum 
entropy principle. Every observable that one can mea- 
sure is interpreted as a constraint. Then one introduces 
Lagrange multipliers and maximizes the entropy (or neg- 
ative information) under these constraints. Using the 
laws of thermodynamics, one shows that the Lagrange 
multipliers are related to the thermodynamic control pa- 
rameters. When one uses only the constraint that the 
average energy of the system has to take on a certain 
value, one ends up with the Boltzmann-Gibbs canonical 
distribution. 

We will now extend these considerations and use the 
maximum entropy principle to obtain the least biased 
distribution for f((3; Xi). As a likelihood estimator we use 
the Shannon entropy, though in principle other choices 
such as the Tsallis entropy [39| are possible as well. The 
entropy associated with the distribution / is 

S(X t ) = -(In/GSjAi)^. (6) 

Clearly the distribution f((3; Xi) has to be normalized. So 
a first property of the distribution /(/?; Xi) that one has 
take into account is (1)^ = 1. Given some complex sys- 
tem in a stationary nonequilibrium state one may have 
additional information on the system which imposes some 
additional constraints. To obtain appropriate constraints 
for superstatistical systems, we briefly repeat the general 
idea of this theory. In each cell, the value of the temper- 
ature is fixed. For the entire nonequilibrium system this 
condition is relaxed and the temperature is allowed to 
vary between the different cells. The crucial assumption 
of superstatistics is that these temperature fluctuations 
have a slow time scale compared with the time scale of re- 
laxation to local equilibrium. The slow fluctuations of the 
temperature cause extra (slow) fluctuations of variables 



3 



like the entropy and the energy in each cell. So it is rea- 
sonable to constrain that the averages of these variables 
should take on certain values. One can still add further 
constraints in terms of some function g(/3), whose pre- 
cise form depends on the nature of the complex system 
considered, i.e. its dynamics, symmetries, and boundary 
conditions. Thus, in the most general case the quantity 
to be optimized is 



S(Xi)-^r(S(0))f>-y(0mb 
-A 3 (5(/3))/3-A 4 (l)^ 



(7) 



with V being an arbitrary constant (taking out a common 
factor out of the definition of Ai and A2 will turn out to 
be useful in the following) . Using the well-known formula 
S(p) = lnZ(/3) + f3E(/3) and renaming (A x + A 2 ) -> A 2 
one obtains 



S{\)-^Z(fi))p-y(fiE{p)h 
-h(g(P))p ~ A 4 (l) /3 . 



(8) 



The optimization of this expression results in the follow- 
ing distribution 



/(/3;Ai) = 



Z(J3) 



-Xi/V 



Z{Xi) 



■ exp 



' V 



(9) 



with Z(Xi) a normalization constant that is fixed by the 
condition = 1. 

We now relate our general result to previous work 
obtained in the literature. In [35| . the authors maximize 
the sum of S(Aj) and (S(0)}p under the constraint of the 
normalization of /(/3;Aj) only. This coincides with our 
approach in case the Lagrange multipliers of expression 
© are chosen in the following way: Ai/V = X2/V = — 1 
and A3 = 0. This results in a distribution that is usu- 
ally not normalizable. For this reason in j3o| the do- 
main of [3 is restricted to a finite range when simple ex- 
amples are studied, such as n non-interacting classical 
Brownian particles. Closely related is also the research 
of Crooks [361 ]. He studies general nonequilibrium sys- 
tems, without assuming that the system can be divided 
into different cells that reach local equilibrium. Crooks 
advocates that instead of trying to obtain the probabil- 
ity distribution of the entire nonequilibrium system one 
has to try to estimate the 'metaprobability', the proba- 
bility of the microstate probability distribution. Crooks 
also uses a maximum entropy principle but puts A3 — 0. 
A main difference is that Crooks does not assume lo- 
cal equilibrium in the cells, hence his approach, though 
an interesting theoretical construction, does not give a 
straightforward physical interpretation to the fluctuat- 
ing parameter (3. The advantage of our approach is that 
one obtains a local fluctuating temperature that coin- 
cides with the thermodynamic temperature and that can 
in principle be measured. The work of Crooks is used by 
Naudts [371 ] to describe equilibrium systems. The author 



shows that some well-known results of equilibrium sta- 
tistical mechanics can be reformulated in a very general 
context with the use of the concepts introduced in [l|, [3(| . 



IV. PHYSICALLY RELEVANT 
DISTRIBUTIONS 

We now discuss some physical properties that should 
be satisfied by the distribution coming out of the entropy 
maximization procedure. Physically one would expect 
the superstatistical distribution /(/?; Ai) to vanish at very 
low and very high temperatures. Assume for the moment 
that no additional constraint exists, i.e. g((3) = 0. In this 
case one can immediately obtain the sign of the various 
Lagrange multipliers by studying the limiting behavior of 
/(/?; Ai) for (3 — > and /3 — ► 00. In the high temperature 
limit, the distribution is proportional to 



Z(/3) 



-Ai/V 



Z(X t 



(10) 



The partition function Z(Ji) usually diverges at high 
temperatures (the entropy becomes infinite). As a con- 
sequence, for physical reasons, the sign of Ai must be 
positive. In the low temperature limit, the energy and 
the entropy go to a constant, lim^oo S(/3) = So and 
lim^_»oo E(f3) — Eo. In this limit, the distribution is 
proportional to 



1 



Z{\) 



exp 



-Ax 



So 

y 



/3(A 2 -A X ) 



Eo 

y 



(ii) 



Therefore, the sign of (A2 — Ai) E /V must be positive. 
Clearly, when a non-trivial function g(0) ^ is imple- 
mented, one has to take into account the limiting be- 
havior of this function as well. For a lot of models 
Eq = 0. In these cases the temperature dependence of 
lim l g_, 00 f(f3; Ai) is solely determined by g{0). This shows 
that implementing a non-trivial function g{(3) 7^ as an 
extra constraint can be important. 

Our reasoning assumes that the low temperature limits 
of S((3) and E((3) are finite constants. This is generally 
true, and is known as the third law of thermodynam- 
ics, but this limit is only taken care of in an appropriate 
way if one uses quantum statistical mechanics. For ex- 
ample, it is well known that the entropy of the classical 
ideal gas diverges at low temperatures. Therefore we will 
now illustrate the general theory with two examples, the 
quantum harmonic oscillator and the classical ideal gas. 
We will come back to the issue of the low temperature 
limit when we study the classical ideal gas. 



SUPERSTATISTICAL QUANTUM 
HARMONIC OSCILLATOR 



As a first example we study n 1-dimensional non- 
interacting quantum harmonic oscillators with tempera- 
ture fluctuations. The Hamiltonian of a single oscillator 
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with mass m and frequency u> is 



H = 



1 

2m 



P 



1 2 2 

— muj x , 
2 



(12) 



with p the momentum operator and x the position oper- 
ator. The energy levels of the oscillator are well-known 
to be 



Ei 



(13) 



with i = 0, 1, 2, . . .. The partition function and the en- 
ergy of the n oscillators become 



Z(J3) = (e h ^/ 2 



E{(3) 



nhjjj 



1 



(14) 



Inserting these formulas into the expression for the dis- 
tribution of the inverse temperature results in 



Z(\i) (1 



exp 



(15) 



with A3 = and V — n. The high and low temperature 
behavior of this distribution is 



Jim f(P;Xi) 



{hwl3) Xl 



(16) 



Clearly, the distribution f(j3; Ai) vanishes at high and 
low temperatures when A2 > Ai > 0. For this quantum 
mechanical example, the low temperature limit of the en- 
ergy Eq = lim^oo E(f3) = hcu/2 is a finite constant. As 
a consequence, no extra constraint (A3 = 0) is necessary 
to obtain a physical relevant distribution. The distribu- 
tion /(/?; Ai) is plotted in Fig. [T]for the example A2 = 2 
and Ai = 1. 



VI. SUPERSTATISTICAL CLASSICAL IDEAL 
GAS 

As a second example we study a 3-dimensional classical 
ideal gas. The gas consists of n particles with mass m 
and is enclosed in a box with a volume equal to unity. 
The partition function and the energy of the ideal gas 
are 



Z((3) 



27T771 



3n/2 



and E(J3) 



3 n 
2j 



(17) 



Inserting these formulas into the expression for the dis- 
tribution J9]) results in 



P'in\ 1 /2V 



exp(-A 3 g(/3)) 



(18) 



0.25 



0.15 




0.05 



FIG. 1: Plot of the distribution of the inverse temperature 
obtained for a set of non-interacting harmonic oscillators. The 
values of the parameters are ?uj = 1, A2 = 2 and Ai = 1. 



The special case X\/V = —1 and A3 = was already 
studied in [35]. As mentioned before, in that case the 
distribution (jT5J) is not normalizable and one has to re- 
strict the values of f3 to a finite range. In [37j the author 
noticed that an inverse gamma distribution is obtained 
for the choice X x /V = -1, g{0) = E{0) and A 3 > 0. 

Let us now comment on physically reasonable choices 
of the function <?(/?). On physical grounds, in the various 
experimental applications of the superstatistics concept 
so far [T^.[l^. [26l]30l [33l.[33|. essentially three relevant dis- 
tributions f((3;Xi) were observed for examples described 
by the superstatistical classical ideal gas: The gamma 
distribution, the inverse gamma distribution and the log- 
normal distribution. Some theoretical reasoning can be 
given why this is so and why the above three distri- 
butions span up three relevant universality classes. It is 
now interesting to see that our generalized maximum en- 
tropy principle, in contrast to previous theoretical work, 
contains all these physically relevant cases. Depending 
on the choice of the function g(/3) and the values of the 
Langrange multipliers A^ one can extract the three rele- 
vant universal distributions out off expression (|18p . For 
convenience, we put V — in/ 2. The gamma distribution 
is obtained for g((3) = (3, Ai > and A 3 > 0: 



/(A Ai) 



Z(Xi) 



exp(-/3|A 3 |). 



(19) 



The inverse gamma distribution is obtained for g{(3) 
Ai < and A 3 > 0: 



mXl) = Z(A0 6XP 







(20) 



The lognormal distribution is obtained for g(f3) = (ln/3) 2 
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and A3 > 0: 



Z{K) 

1 1 

Z'(Xi) ~p 



exp 1 



|A 3 |(ln/3) 2 



exp 



-|A3|(ln/3-A 4 ) 2 ), (21) 



10" 



with 

A 4 = ^(A 1 + 1). (22) 

Unlike the quantum mechanical case, for classical com- 
plex systems usually g{(3) 7^ is needed to make expecta- 
tions formed with /(/3; Xi) converge. This function <?(/3) 
is determined by additional information that one has on 
the complex system under consideration (an example will 
be treated in the next section). 

Unlike the quantum mechanical case, for the classical 
ideal gas one has to be careful with a range of f3 that goes 
from to 00. For this example the limiting behavior of 
the energy and the entropy at low temperatures is 

lim E(0) = 0, and lim S(j3) = -00. (23) 

(3 — >oc /5 — >oo 

The limit of the energy is acceptable from a thermody- 
namical point of view, the limit of the entropy is not. 
Clearly, the problem arises from the fact that the clas- 
sical treatment of an ideal gas in the low-temperature 
limit does not make sense, one certainly has to take quan- 
tum corrections into account. However, when /(/3; Xi) is 
vanishing in this limit, the contribution of the quantum 
region to the average values of the observables will be 
negligible. Notice that the three aforementioned distri- 
butions (fT§|) . (|2"U)) . (f2"Tj) all have a single peak at a well 
defined temperature. So as long as this single peak is 
situated in the classical region one can use classical mod- 
els in the context of superstatistics although one has to 
be careful in evaluating the low temperature behavior of 
/(# Xi) itself. 

VII. TURBULENT TAYLOR-COUETTE FLOW 

As a final example we now apply our methods to a com- 
plex system that is not analytically solvable anymore: 
Turbulent Taylor-Couette flow. Ideally, for a supersta- 
tistical statistical mechanics description of this system, 
one would measure the set of all positions and velocities 
of a large number of test particles in the flow. This is 
not possible and hence, as in previous papers [2(], we re- 
strict ourselves to the information that one can get out 
of a scalar time series, a single measured velocity compo- 
nent v(t) as a function of time t. We use data from an 
experiment performed by Lewis and Swinney [40| . The 
stationary probability distribution p(u) of the velocity 
difference u(t) = v(t + 8) — v(t) at a given scale 8 is well- 
known to exhibit non-Gaussian behavior, see Fig. [5] for 
an example. 
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FIG. 2: Stationary distribution p(u) of velocity differences 
u(t) as measured in the experiment of Swinney et al. [ipl ] at 
Reynolds number Re = 69000 and scale 8 — 64. The dashed- 
dotted line is a Gaussian distribution 0.1280 exp(—0.0515it 2 ), 
whereas the solid line corresponds to the superstatistical for- 
mula (f24)l . 




FIG. 3: Example of a probability distribution /(/?; Ai) as ex- 
tracted from the measured turbulent time series of velocity 
differences for Re = 69000 and 5 — 64. The solid line is 
a fit to the lognormal distribution (|21|l . with A3 = 3.8516, 
A 4 = -2.303 and Z'(Ai) = ^n/X 3 . 



It has been previously shown that superstatistical tech- 
niques can be successfully used to model the statistics 
of turbulent velocity fluctuations 0, [2(J HH, [26[ . For a 
measured time series u(t), the parameter (3 simply corre- 
sponds to a local inverse variance of the measured signal, 
and the 'cells' of the superstatistics approach correspond 
to time slices of a suitable length where this variance 
is measured. The turbulent velocities are well approxi- 
mated by the model of a classical ideal superstatistical 
gas, meaning that for certain time intervals the signal 
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FIG. 4: Dependence of the Lagrange multipliers A3 and A4 
on the scale 8 and Reynolds number Re, Re = 540000 (solid 
lines), Re = 266000 (dotted lines), Re = 133000 (dashed- 
dotted lines) and Re = 69000 (dashed lines). The inset shows 
a magnification for large values of 5. 



is Gaussian with a given variance, then it changes to 
another Gaussian with a different variance, and so on. 
The validity of the above approximation and the neces- 
sary time scale separation has been checked in a previous 
paper In that paper also a general method was in- 
troduced how to to extract the relevant time slicing (the 
superstatistical cell size) and how to extract the distri- 
butions Aj) from the signal. We do not describe this 
here in detail, but refer to the paper Using these 
techniques, we determined the distribution /(/3; \) from 
the measured time series, using the experimental data of 
Swinney et al. for various scales S and Reynolds numbers 
Rc. In all lognormal distribution turns out to be 

a reasonable fit for the experimentally observed distribu- 
tion f(f3;\i), see Fig. [3] for an example. However, the 
parameters of this lognormal distribution depend on 5 
and Rc in a nontrivial way. Our results are summarized 
in Fig. |H 



The relevance of lognormal distributions is to be ex- 
pected due to the multiplicative random processes un- 
derlying the fluctuating energy dissipation in turbulent 
flows. In other words, the cascade picture of turbu- 
lence suggests that the constraint g{0) in the general- 
ized maximum entropy principle should be of the form 
g{0) = (ln/3) 2 , leading to lo gnormal distributions. More 
surprising is the fact that our data analysis indicates that 
there is a distinguished scale In 5* « 3.5, or 6* » 32, 
where the obtained fitting parameters A3,A4 are inde- 
pendent of Reynolds number. For 6 < 5*, A3 increases 
with increasing Reynolds number, whereas for 5 > 6* it 
decreases. A4 shows the opposite behavior, it decreases 
with Re for 5 < 5* and increases for 5 > S* . 

One may check the quality of the superstatistical 
model approximation by numerically evaluating the dis- 
tribution Q 

p(u)« r m\i)\ffe-iP u2 dl3 (24) 
Jo V 2vr 

and comparing it with the measured stationary distribu- 
tion p(u). Here /(/3; Ai) is a lognormal distribution with 
parameters as given in Fig. 2] The solid line in Fig. [5] 
shows this curve (f24f for the example Re = 69000, 5 = 64. 
Clearly, there is an excellent agreement between the ex- 
perimentally measured distribution p(u) and the super- 
statistical approximation (p?4"|) . 

Our turbulence example illustrates that the Lagrange 
multipliers in the generalized entropy maximization prin- 
ciple, A3 and A4, do have physical meaning. Under dif- 
ferent conditions, in our case fixed by the scale 5 un- 
der consideration as well as the Reynolds number of the 
flow, these intensive parameters have different values (see 
Fig. 0J. In fact, one could go as far as regarding the re- 
sults in Fig. [4] to describe a kind of 'equation of state' 
of the turbulent Taylor-Couette flow, providing the de- 
pendence of the intensive parameters A3 and A4 on given 
parameters of the flow pattern, such as scale 8 and Rc. 
All this illustrates that the generalized maximum entropy 
principles developed in this paper are not only a mathe- 
matical exercise, but of true physical relevance for a va- 
riety of classes of complex systems, when only a certain 
limited amount of information on the system is available. 

VIII. CONCLUSION 

In this paper we developed a maximum entropy prin- 
ciple for superstatistical systems of various kinds. This 
principle allows for the determination of the superstatis- 
tical distribution /(/3; Aj) of the fluctuating intensive pa- 
rameter 0, given some prior information on the complex 
system under consideration. Our formalism further de- 
velops previous work of Abe et al., Crooks, and Naudts, 
and contains physically relevant superstatistical univer- 
sality classes, such as lognormal superstatistics, gamma 
superstatistics and inverse gamma superstatistics, as spe- 
cial cases. We dealt with 3 important physical examples, 
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the superstatistical quantum harmonic oscillator, the su- 
perstatistical classical ideal gas and time series as gener- 
ated by a turbulent Taylor- Couette flow. For the quan- 
tum case, a new single-peaked distribution Xi) as 
displayed in Fig. [T] arises quite naturally out of our max- 
imum entropy approach, whose physical relevance can 



be checked in future experiments. For classical systems, 
other types of distributions are relevant, such as the log- 
normal distribution for turbulent flows, as displayed in 
Fig. EI Our approach is a further step to arrive at a gen- 
eralized statistical formalism relevant for large classes of 
complex systems with time scale separation. 
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